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1 .  INTRODUCTION 


The  present  University  of  Michigan  surface  field  measurement  facility 
[1]  covers  the  frequency  range  from  100  MHz  to  4770  MHz,  corresponding  to  a 
wavelength  of  3  m  at  the  lowest  frequency  and  6.4  cm  at  the 
highest  frequency.  However,  at  frequencies  below  120  MHz  the  performance 
of  the  chamber  and  its  associated  equipment  deteriorates  rapidly.  For 
accurate  measurements  it  is  necessary  that  the  sensor  be  small  compared 
with  the  relevant  dimensions  of  the  target,  and  since  there  is  a  minimum 
size  of  sensor  that  can  be  constructed  and  that  has  sufficient  sensitivity, 
there  is  in  turn  a  minimum  size  target  that  can  be  used.  In  practice, 
therefore,  the  target  must  have  overall  dimensions  greater  than  (about) 

25  cm,  implying  a  resonant  wavelength  of  0.5  m  (f  =  600  MHz).  The  net 
result  is  that  the  present  facility  enables  us  to  measure  down  to,  at 
best,  a  frequency  of  0.2  f^,  where  f^  is  the  lowest  resonant  frequency 
of  the  target.  This  is  not  quite  low  enough  to  define  the  low  frequency 
behavior  of  the  target.  It  is  believed  that  a  frequency  of  0.1  f^  or, 
preferably,  0.05  f^  must  be  attained  to  specify  the  low  frequency 
asymptote,  and  since  it  is  inconceivable  that  any  simple  modification  of 
the  present  facility  would  enable  us  to  get  down  to  30  MHz,  the 
construction  of  a  separate  facility  for  measuring  the  surface  fields  at 
one  (or  more)  frequencies  f  _<  0.05  f^  was  undertaken. 

The  main  purpose  of  the  new  facility  is  to  extend  the  presently 
measured  curves  of  surface  magnetic  field  (or  current)  down  to  frequencies 
of  30  MHz  or  below  by  specifying  the  low  frequency  asymptote  in  each 
case.  A  measurement  at  a  single  frequency  would  suffice  and,  in 


principle  at  least,  the  lower  tne  frequency  tne  better.  Since  the  m^:gnetic 
field  is  of  interest,  it  is  natural  to  consider  a  quasi -magnetostatic 
facility  consisting,  for  example,  of  a  pair  of  Helmholtz  coils,  and  it 
is  desirable  that  the  facility  be  able  to  use  models  similar  (or,  indeed, 
identical)  to  the  ones  employed  in  the  present  facility. 

Theoretical  analyses  of  the  fields  produced  by  dc  and  rf  excited 
Helmholtz  coi's  are  presented  in  Section  2,  and  other  pertinent  studies  such 
as  the  effect  of  a  non-perfectly  conducting  or  non-magnetic  test  body  are 
also  addressed  here.  Section  3  describes  the  design  and  implementation 
of  the  two  facilities  that  were  evaluated--one  of  plywood  construction 
and  the  other  all  metal.  Field  maps  and  the  frequency  response,  are 
presented  in  Section  4  and  compared  with  the  design  criteria.  Section  5 
gives  sample  measurements  for  a  sphere,  cylinders  and  an  F-106B  model 


aircraft. 


2.  THEORETICAL  CONSIDERATIONS 


2.1  Analysis  of  Fields  Produced  by  dc  Excited  Helmholtz  Coils 

2.1.1  Introduction.  The  Helmholtz  coil  arrangement  consists  of  two 
circular  coaxial  coils,  each  of  the  same  mean  radius  a,  spaced  apart 
(between  the  midplanes)  by  the  distance  a,  and  with  conductor  (or  wire) 
diameter  kept  to  a  suitably  small  fraction  of  a  [2].  When  the  coils  are 
excited  by  direct  currents  in  additive  series,  i.e.,  when  each  coil  carries 
a  similarly-oriented  steady  current  I,  a  near-uniform  axially  directed 
magnetic  field  is  produced  in  the  central  section  of  the  axis  of  the 
system  [3].  In  this  Section  we  develop  the  exact  expressions  for  the 
various  components  of  the  field  produced  by  the  dc-excited  Helmholtz 
coils,  and  derive  some  approximate  expressions  to  examine  the  uniformity 
of  the  fields  in  the  central  region  of  the  system. 


2.1.2  Field  Produced  by  a  Circular  Current.  Let  us  consider  a 
circular  loop  of  radius  a  carrying  a  steady  current  I  with  the  plane 
of  the  loop  oriented  in  the  x-y  plane  such  that  its  center  coincides 
with  the  origin  of  a  cylindrical  coordinate  system  (r,  (j),  z).  Since  the 
exact  expressions  for  the  field  components  produced  by  such  a  system  are 
well-known,  we  shall  simply  quote  the  relevant  expressions  (see  [4,5] 
for  the  detailed  derivations). 

The  magnetic  vector  potential  at  a  field  point  P(r,z)  is 
entirely  4>-clirected  and  is  given  by 
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A.  •  #  <" - ^  ■  [2.U 

’■  ^  '  (a-  +  r-  +  z-^  -  2ar  cos  -i) 

where  u  is  the  permeability  of  the  medium  and  the  other  quantities 
are  as  defined  earlier.  Equation  (2.1)  can  be  expressed  in  terms 
of  known  integrals  as 

A  =  lii  (  i  )’'-[(!  -  i<‘)  K  -  E]  (2.2) 

■P  ir<  ' 


where 


4ar _ 

C(a  +  r)2  +  z2] 


(2.3) 


and  K  and  E  are  complete  elliptic  integrals  of  the  first  and 
second  kinds,  defined  as  [6]; 


n/2 

K  =  f  - - - ^  . 

^  (1  sin2  6)'/' 

^/2 

E  =  f  (1  -  K  ^  sin^e )  ^  de  . 


(2.4) 


(2.5) 


Assuming  u  to  be  constant,  the  two  field  components  can  now  be 
obtained  from 


H  = 
r 


1 


u 


3Z 


H 

z 


ur 


(2.6) 


Us  i  ng 


ihroLiyh  (2.6)  the  following  two  exact  expressions  are 


obtained  for  the  field  components  [3,4]; 


V  ■  27 


rC(a  +  r)2  + 


[(a  +  r)2  +  z2]’ 


+  a2  +  ,’2  +  2^  ^ 

{a  -  r)2  +  z2  . 


K  ,  -  .rz. 

(a  -  r)2  +  z2 


(2.7) 


(2.8) 


It  can  be  shown  (see  Appendix  A)  that  near  the  axis  (i.e.,  as 
r  -►  0)  the  fields  may  be  approximated  as 


^  C(a  +  r)2  +  z2]^/^ 


(2.9) 


H  =  ^  1  +  I  ai  -,..4zi_  ^2  +  o(r3)  .  (2.10) 

"  2  (^2  +  22)3/2  [_  4  (a2  +  z2)  J 

From  (2.9)  and  (2.10)  it  is  found  that  the  fields  on  the  axis  (r  =  0) 


=  0  , 


_  la^ 


f“  - ^ - 'o  3/2  ~  ^ 

^  (a"  +  z^y'^  ^ 


(2.11) 


and  as  shown  in  Appendix  A,  these  are  exact.  -  The  first  three 
derivatives  with  respect  to  z  of  the  function  f(z)  defined  in  (2.11)  are 


f'(z)  =  - 


f"(2) 


(a^  +  z2)3/2 

.  3(a2  -  4z^) 

’  (a2  +  z2)2/2 


(2.12) 


Using  (2.11)  and  (2.12)  it  is  seen  that  with  two  similar  coils 
located  at  z  =  ta/2  the  first  three  derivatives  of  the  total  field 
(H^  in  this  case)  on  the  axis  vanish  at  z  =  0.  This  is  the  basis 
for  the  design  of  Helmholtz  coils  discussed  in  the  next  section. 

2.1.3  Helmholtz  Coils. 

2. 1.3.1  Fields  on  the  Axis  (r  =  0)  of  Two  Coils.  For 
generality  we  consider  two  circular  coils,  located  at  z  =  ±h  (referred 
to  as  coils  1  and  2,  respectively),  each  of  radius  a  and  carrying  a 
steady  current  I.  Using  (2.11)  it  can  be  shown  that  the  exact  field 
for  r  =  0  produced  by  the  system  are 

(2.13) 

Ja2  +  (z-h)2]^/^  [a2  +  (z+h)2]3/2_ 


Differentiating  H^  with  respect  to  z,  the  first  four  derivatives  at 
z  =  0  are  found  to  be 


0 


a^H. 


3Z' 


=  -3Ia' 


4h^ 


(a2  +  h2)^/2 


(2.14) 
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wv 


h**-*...  .t—  ^kim  A»«  -1^  %*»  >♦«  !.»•  tw  ■fcm  ac 


3Ia^ 

(a2  +  h2)5'2 


•ISa**  ■^  180a%^  -  120h‘*) 
(a2  +  h2)3 


It  is  clear  from  (2.14)  that  for  h  =  a/2  the  first  three  derivatives 
of  vanish  at  z  =  0  and 


^“*^2  u  144  4! 

3Z-  ■  “  '25  • 


at  z  =  0  , 


(2.15) 


where 


»o  ■  "z 


at  2=0  5^'^a 


(2.16) 


is  the  value  of  at  the  center  of  the  Helmholtz  coils.  For  later 
use  we  now  write  the  following  Taylor  series  approximation  to 
near  z  =  0  obtained  by  using  (2.13)  through  (2.16): 


(z/a)‘ 


(2.17) 


for  -  1/2  <  z/a  <  1/2. 

Equation  (2.17)  is  convenient  for  estimating  the  variation  of 
the  field  on  the  axis  near  the  central  region  of  the  Helmholtz 
coils.  From  now  on  we  shall  assume  that  h  =  a/2,  i.e.,  the 
separation  distance  between  the  two  coils  is  equal  to  a. 


2. 1.3.2  Fields  in  the  Central  Region  (r  0).  Using  (2.9) 
and  (2.10)  the  fields  in  the  central  region  of  the  Helmholtz  coils  may 
be  approximated  as 


1  2 

"z  '  «z  *"z 

1  2 


(2.18) 

(2.19) 


with 


3Ia2  (z  +  a/2)r _ 

C(a  +  r)2  +  (2  ;  a/2)2]^/2 


(2.20) 


H 


z 


1 

2 


Ia2  _ ] _ 

[a2  +  (z  +  a/2)2]3'2 


1  +  3  a2  -  4(z  ;  a/2)2  ^2  , 

[a2  +  (z  ;  a/2)2]2  J 


(2.21) 


where  the  subscripts  1  and  2  on  the  left-hand  sides  of  (2.20)  and  (2.21) 
correspond  to  the  neQative  and  positive  signs,  respectively,  on  the 
right-hand  sides.  From  (2.18)  and  (2.21)  it  can  be  seen  that  in  the 
z  =  0  plane. 


where  is  given  by  (2.16).  Equation  (2.22)  shows  that  the  Helmholtz 
coils  maintain  a  constant  and  axially  directed  field  in  a  plane 
parallel  to  the  coils  and  passing  through  the  center. 
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We  are  now  in  a  position  to  obtain  some  approximate  expressions 
to  estimate  the  uniformity  of  the  central  region  fields.  Examination 
of  (2.18)  through  (2.22)  reveals  that  exactly  at  the  center  and  in  the 
plane  z  =  Q,  =  0.  On  the  axis  (r  =  0),  =  0  and  the 

variation  of  may  be  estimated  by  using  (2.17)  which  is  rewritten 


Hz  =  H,  1 


(2.23) 


where 


1  z  1 

2  a  2 


(2.24) 


If  desired,  the  variation  of  with  r  in  this  region  may  be  found 

using  (2.21).  From  (2.18)  and  (2.20)  the  radial  component  of  the  field 
is  given  by 


1  \ 


'  *  *  [I  -  1)^1 


1  +  -’> 

>  a/ 


.rf  .  '  ^  s/: 
^  J  • 


(2.25) 


1  z  1 

■  2  a  2 


and  from  (2.24)  and  (2.25)  it  can  be  shown  that 


0.5  X  10 


t  ^  and  r 


=  1.0  X  10' 

Ho 


at  -  = 
a 


.o 
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More  accurate  and  detailed  expressions  for  the  variation  of  and 
inside  the  coils  can  be  found  using  the  exact  formulas  for  the  fields 
given  in  the  next  section. 


2. 1.3.3  Exact  Expressions  for  the  Fields.  Using  (2.7)  and 
(2.8)  the  following  exact  expressions  for  the  field  produced  by  the 
Helmholtz  coils  can  be  obtained: 


(2.26) 


(2.27) 


V  2t 
n 


r[(a  *  r)2  . 


a2  +  +  z2 

,K  +  - 0-  E  ,  (2.28) 

(a  -  r)2  +  z-  J 


(2  -  r2  -  72 


- ! -  K  +  - —  E 

[(a  ♦  r)^  ♦  zn’"  .  "  (a  -  r)^  +  z2  " 


(2.29! 


where  z  =  z  -  a/2,  z  =  z  +  a/2, 

K  ,K  are  the  elliptic  integrals  defined  by  (2.4)  with  arguments 

■:  ,<  obtained  from  (2.3)  after  replacing  z  by  z  -  a/2, 

1  2 

z  +•  a/2,  respectively;  and 

E  ,E  are  the  elliptic  integrals  defined  by  (2.5)  with  arguments 

1  2 

■-  obtained  in  a  similar  manner. 

1  2 

It  should  be  noted  that  for  r  =  0.  H  =0  for  all  z. 


2.1.4  Numerical  Computations.  Equations  (2.26)  and  (2.27)  were 
programmed  to  compute  the  radial  and  axial  components  of  the  fields 
normalized  with  respect  to  the  field  at  the  center,  i.e.,  H^/H^  and 
respectively,  as  functions  of  the  normalized  coordinates 
r/a  and  z/a,  where  a  is  the  radius  of  each  coil.  Selected  results  for 
H^/H^  and  H^/H^  at  r/a  =  0(0.25)0.75  and  z/a  =  0(0.25)0.5  are  given  in 
Table  2.1.  The  computer  program  is  listed  in  Appendix  A. 

The  computations  were  carried  out  for  coils  separated  in  distance 
a  where  a  is  the  radius  of  the  coil.  As  expected  and  seen  from  the 
Table,  the  component  is  dominant.  At  the  center  (r/a  =  0,  z/a  =  0) 


Table  2.1 

Numerical  Values  for  Normalized  Field  Components 


r/a 

Z/a 

Hr/% 

0 . 000000 

0.0 

0.0 

1 .0000000000 

0.000000 

0 . 250000 

0.0 

0.9958016347 

0.000000 

0.500000 

0.0 

0.9458241852 

0 . 2 S 0000 

0 . 0 

0.0 

0.9982148100 

0. 2fi0000 

0 .  :'50ooo 

0.00305210311 

1.0069158788 

.  2ti0C00 

0.500000 

0.045B2 72520 

0.9718178124 

w .  4V>yV99 

O.u 

0.0 

0.9666290499 

0 . 4yV999 

0 . 250000 

-0.037*608067 

1 .0361427785 

0 . 4VVV99 

0.500000 

0.0877231827 

1.0812025957 

0  .  .M999V 

0 . 0 

0.0 

0.7981811608 

0  .  /-(V999 

0.250000 

-0.2^57322.9654 

1 .0205464277 

'  .  ,M9V99 

0.500000 

0. J 179854157 

1.5052027031 

it  is  unity  along  the  z-axis  and  decreases  to  0.945  (five  percent)  at 
z/a  =  0.5  (r/a  =  0)  and  in  the  horizontal  plane  through  z/a  =  0  decreases 
to  0.967  (four  percent)  at  r/a  =  0.5.  Note  at  r/a  =  0.5,  z/a  =  0.5  the 
field  increases  to  1.081  which  is  expected  since  this  point  is  closer  to 
the  conductor  (coil)  than  the  other  two  points.  Observing  also  that  at 
r/a  =  0.25,  z/a  =0.25  the  field  is  1.007,  it  is  reasonable  to  conclude 
that  for  a  spherical  volume  of  r/a  =  0.5  located  at  the  center  of  the 
coil  system  the  field  will  be  within  ±5  percent  of  the  center  value  H^. 

On  the  other  hand,  in  a  smaller  circular  region  defined  by 

z/a  =  ±0.25  ,  r/a  =  0.25 


the  deviation  of  the  field  is  less  than  one  percent  from  that  at  the  center. 
For  the  dimensions  selected  (see  Section  3)  this  corresponds  to  a  cylindrical 
volume  30  cm  in  diameter  and  30  cm  in  height.  The  scale  models  typically 
used  in  our  measurements  are  no  larger  than  this. 

2 . 2  Analysis  of  Fields  Produced  by  ac  Excited  Helmholtz  Coils 

The  standard  Helmholtz  coil  arrangement  consists  of  two  circular 
coaxial  coils  of  radius  a  spaced  a  apart.  If  each  coil  carries  a 
similarly  directed  constant  current  I  then,  as  shown  in  Section  2.1.4, 
there  is  a  cylindrical  region 


where  the  radial  magnetic  field  is  negligible  and  the  longitudinal 
magnetic  field  is  constant  to  within  one  percent.  It  is  proposed  that 


this  constitutes  the  test  region  of  our  planned  low  frequency  (magnetic) 
facility.  In  order  to  accommodate  a  model  of  overall  dimensions  up  to 
0.25  m,  it  is  necessary  that  the  radius  of  the  coils  be  at  least  0.5  m, 
and  the  proposal  is  that  a  =  0.61  m  (2  feet). 

There  are  two  sources  of  a  frequency-dependent  non-uniformity: 

(i)  variations  in  the  amplitude  and/or  phase  of  the  current  I  excited 
in  each  coil,  and  (ii)  the  effect  of  the  phasor  addition  of  the 
contribution  of  even  a  uniform  current  on  the  near  (and  far)  field  of 
a  coil.  Since  (i)  is  a  function  of  the  manner  in  which  the  coils  are 
excited,  we  shall  concentrate  only  on  (ii). 

Following  Section  2.1.2,  consider  a  circular  coil  of  radius  a  in 
the  plane  z  =  0  of  a  cylindrical  polar  coordinate  system  (r,(j),2).  If 
the  coil  carries  a  constant  (uniform)  current  I,  the  vector  potential 
at  an  arbitrary  field  point  is  A  =  with 

<p 


ik  /a2+r2+z2-2ar  cos  a 


1  -  Hli 

\  ~  2tt 


cos  a  e 


/a2+r2+z^-2ar  cos  a 


(2.30) 


where  a  time  dependence  exp(-i(jjt)  has  been  assumed  and  suppressed. 
For  small  ka,  kr,  kz,  the  exponential  can  be  expanded  to  give 


cos  a/a2+r2+z2-2ar  cos  a  da  + - 


(2.31) 


implying 


A  = 


3<2 


1(1  -  k2)K  .  1  .  ^  E 


Since ,  for  <2  <<  1 ^ 


T  1  +  0  {<'*) 


E  = 


J  jl  -  ^  <2  +  Q  (^4)J 


it  follows  that  A  =  0  when  <  =0,  as  required. 
From  (2.33)  and  (2.35) 


r  cos  a/a2+  r2+22  -  2ar  cos  a  da  = 


®  /ar  j(l  -<2)k  - 


3<2 


I 


(2. 


(2 


2 


(2 


and  since 


f 


cos  a  da  =  0  , 


we  have 


=  xJT  (’-r) 


K  -  E  - 


2ar  /k 


3  V< 


!  ( 1  -  <2 


-  n  -  piE^  +  0  (k3) 


.  (2 


The  first  term,  independent  of  k,  is  identical  to  that  in  Eq.  (2.2) 


In  terms  of  A  the  radiated  field  is 


H  =  j  V  X  A  ,  E 


i  k  Z  1  A 
y 


where  Z  =  1/Y  is  the  free  space  impedance,  implying 


"r  =  -  ^  ^  •  "z  =  If  • 


YE^  =  ik  -  A,  . 


(2.39) 


with  all  other  components  zero.  It  follows  immediately  from  (2.38)  that 


YE 


a2  +  r^  +  z^ 


ik  j"(a  +  r)2  + 

2*'^  [  J  L(a  *  >-)2  +  z= 


K  -  E 


k2 

6 


r. 


-  ^  < 


1 


(a  -  r)2  +  z2k  -  (a2  +  r2  +  z2)E  +0(k3)j  (2. AO) 


and  using  the  expansions  (2.36)  it  can  be  shown  that  =  0  on  the  axis 
r  =  0  of  the  coi 1 . 

To  determine  the  magnetic  field  we  first  note  that 


If.  ’E—  a2  -  r2  +  z2  ^ 

■  2'"  (a  a  rF  *  z^  ’ 


<^Z 

Aar 


Also,  from  [7] 


implying 


^  ,  i 

die  ,  < 

L 


I' 


<  ■'  1/ 

rj 


- 1-  Jk  . 

-2  j  1  .  .2 


L 


'^1  -  5-)  K  -  E) 


After  some  rather  tedious  manipulation  it  now  follows  that 


and 


I 

Zirr 


{ (a+r)2+z2} 


ITT 


a2  +  p2+22  1^2  r 

+  £  (a2+r2+z2)K 

(a-r)2+z2  2  ^ 

-  ((a-r)2+22]E]  +  0  (k^)  (2.41) 


=  i- 


2  2it  {(a+r)2+z2}i/2  |_  (a.r)2+z2 


K  +  aiiliilL  E  +  ^  ]  (a2.r2.22)K 


+  ^{a+r)2+z2j  E  I  +  0  (k3) 


(2.42) 


The  leading  terms  in  (2.31)  and  (2.32)  are  in  agreement  with  those  given 
in  Eqs.  (2.7)  and  (2.8). 

For  small  r  such  that  2ar  <<  a2+z2  we  have,  to  the  first 
non-zero  term  in  r. 


k2 

2  (a2+z2)  +  0 (k3) 


IL  a^rz 

2 


k2 

1  +  J  (a2+z2)  +  0 (k3) 


(2.43) 

(2.44) 
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1  a- 

2  (a- 


k2 

4-  ^  (a2+z2)  ^  o(k’) 

J 


12, 


E  and  H  are  botfi  proportional  to  r  and  therefore  vanish  on  the 
b  I" 

axis,  whereas  H  is  non-zero  when  r  =  0.  For  H  and  H  the  frequency- 
z  r  z 

dependent  correction  is  proportional  to  k2(a2+z2),  and  in  the 
cylindrical  region  [z]  a  the  correction  is  less  than  one  percent 
if  ka  <  0.1.  This  condition  is  satisfied  for  a  coil  not  exceeding 
0.95  m  in  diameter  operated  at  10  MHz. 

The  proposed  facility  consists  of  two  identical  coils  spaced 
a  apart,  identically  excited.  In  the  resulting  test  region  (of 
overall  length  a/2)  the  variation  of  along  the  axis  can  be 
obtained  from  (15).  The  center  of  the  region  corresponds  to 
z  =  a/2  in  (15),  and  here 


H 


z 


z=a/2 


^  0.3578[1  +  0.625(ka)^  +  0  (k^)] 

Q 


Simi  larly , 


z=a/4 


=  i  0.4565[1  +  0.531(ka)^  +  0  (k^)] 

Q 


and 


z=3a/4 


=  I  0.2560[1  +  0.781(ka)^  +  0(k3)] 

a 


the  sum  of  which  is 


I  0.7125[1  +  0.621(ka)"  +  0  (k^)] 

Q 
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This  differs  by  less  than  one  percent  from  twice  the  value  at  z  =  a/2 
(the  factor  2  corresponding  to  ttie  two  coils).  The  uniformity  of  the 
field  throughout  the  test  region  that  was  derived  in  Section  2.1  for 
static  fields  is  therefore  maintained  if  ka  0.1. 

2.3  Other  Considerations 

In  using  the  magnetostatic  facility  certain  problems  now  arise, 
three  of  which  are  as  follows:  (i)  The  desired  surface  magnetic  field 
is  that  at  a  specific  point  on  the  target  when  illuminated  by  a  plane 
linearly  polarized  wave  incident  in  a  given  direction,  and  this 
excitation  is  quite  different  from  that  provided  by  a  pair  of  Helmholtz 
coils;  (ii)  at  low  frequencies  the  metallized  target  no  longer  looks 
perfectly  conducting;  and  (iii)  the  low  frequency  expansion  is  not 
uniform  in  the  magnetic  properties  of  the  target.  A  non-magnetic  body 
having  u  =  has  no  effect  on  a  magnetostatic  field  (the  surface  magnetic 
field  is  that  of  the  incident  field  alone).  The  static  result  for  a 
"perfect  conductor"  corresponds  to  the  non-physical  situation  for  which 
u  =  0,  and  with  any  actual  highly  conducting  material  there  is  a  frequency 
below  which  the  field  starts  to  diverge  from  that  for  a  perfect  conductor. 
We  shall  examine  these  topics  in  the  reverse  order. 

2.3.1  Quasi-Magnetostatic  Behavior.  To  examine  the  low  frequency 
behavior  of  the  surface  magnetic  field,  consider  the  simple  problem  of  a 
homogeneous  sphere  of  permittivity  e,  permeability  j  and  radius  a 
illuminated  by  the  nlane  electromagnetic  wave 


■  me 


.  i  kz 
X  e 


Z  H 


me 


.  ikz 
y  e 
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Note  that  the  time  convention  (e"^^^)  and  the  direction  of  incidence 


differ  from  those  used  in  Appendix  C,  Part  3. 

As  shown  by,  for  example,  Stratton  [2],  the  total  magnetic 
field  is 


-s 

n=i 


•  n  2n  +  1 
n(n  +  1 ) 


ISp)  *  -  iP)  + 

em  n  ein 


i(n 


(0 

oin 


afip)): 

n  oin' ' 


(2.46) 


where 


a 


n 


and 


b 


n 


u/^^  Va) 

-  r^^^(kia) 

j„(ka) 

u^r^3)(|^a) 

-  r^i^(kja) 

h‘‘'(ka) 

Epf^^^ka) 

-  r^^^(kja) 

jp(ka) 

e^r^^^ka) 

-  r^^^(kja) 

hj,^^(ka) 

(2.47) 


(2.48) 


with  £  =  ^  ~  ^r^o’  and 


r 


(1) 


(x) 


[xj^(x)] 


r^’'(x) 


(2. 


If  'ka  ,  jk^aj  <<  1,  the  coefficients  a^  and  b^  can  be  expanded  in 
powers  of  ka.  In  particular. 


ai 


Up  -  I  5 

-Sr7+0[(ka)  ] 

Up 


implying 
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+  0  [(ka )  ]  , 


implying 


T5  2u  +  3  ^ 


b^  15  2e^  +  3  ^ 


The  low  frequency  expansions  of  the  tangential  components  of 
the  total  magnetic  field  on  the  surface  r  =  a  of  the  sphere  are  then 
as  follows: 


ZHg  =  sin  i>}— V?  ®  T  T~T~2  ^  2e  +  0  [(ka)  ] 


cos  ({!' — ^  cos  8 

Up  +  2 


(2.50) 


+  2  ^  +  0  [(ka)  ] 


We  observe  that  the  static  (zeroth  order)  terms  are  functions  only 
of  the  relative  permeability  and  that  for  a  non-magnetic  material 


(Up  =  1) 


H  = 


independent  of  the  electrical  conductivity.  On  the  other  hand,  for 
a  perfectly  conducting  sphere  it  is  known  that 


ZH, 


2  sin  (})  |cos  0  +  ika  f  20j  +  0  [(ka)  ] 

cl 

t  0  [(ka)  ] 


1 


ZH,  =  T5-  cos  *  <1  +  ika  cos  6 1-1  +  ^ 


[- 


(2.51) 


which  is  otherwise  obtained  by  putting  =  »  and  =  0  in  (2.50) .  We 

remark  that  the  assumptions  made  in  the  expansions  leading  to  (2.50) 

may  be  violated  i^  e  =  “. 

r 

For  a  conducting  material 
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r 


+  iZ 


a 

k 


(2.52) 


where  is  the  (real)  relative  permittivity  and  o  is  the  conductivity. 
The  relevant  factor  is  a/k,  which  is  infinite  for  a  perfect  conductor. 
Moreover,  any  material  for  which  a  ii*  0  "looks"  perfectly  conducting 
in  the  static  limit,  and  if  a/k  is  large. 


k 

1 


(1  +  i 


kZu^a 


(2.53) 


Consider  now  the  exact  expressions  (2.47)  and  (2.48)  for  the 

coefficients  in  the  Mie  series.  We  observe  that  these  involve  k  only 
/  \  1 
via  the  term  r'^'(k  a),  and  for  large  |x| 

1 

r^^^x)  =  1  +  X  =  X  tan  j^x  -  (n  +  1)  fj  ^  °  (|xr^)j  . 


For  large  Im  x  regardless  of  Re  x. 


tan  'x  -  (n  +  l)n/2;  '  i 


implyi ng 


r^‘^(ka)  ika  (2.54 

1  1 

which  tends  to  infinity  as  ku^o  -<■  ®.  Hence,  if  ku^  f  0,  the 
perfectly  conducting  limit  of  (2)  is 


jn(ka) 

h^^^ka) 


(2.55 


as  used  in  the  derivation  of  (6).  We  remark  that  the  same  result 
is  obtained  by  putting  =  0  in  (2). 

For  the  coefficient  b^  the  relevant  quantity  is  e^/r^^^(k  a) 
If  the  body  is  conducting  (a  ^  0),  (7)  and  (9)  imply 


r^^^(k^a) 


,  e 


-iiT/4 


ka 


(2.56 


which  becomes  infinite  as  o  -*■  ®  (perfect  conductivity),  and/or 
k  -  0  (static),  and/or  -*■  0.  In  each  case  (2.48)  reduces  to 


b  >  -  (2.6 

n  - - - 

[kah^^^^ka)]' 

which  is  the  known  result  for  perfect  conductivity,  as  used  in  the 
derivation  of  (2.51) . 

The  conclusions  of  the  above  analysis  are  as  follows.  For 

a  conducting  material  having  a  ^  0,  the  low  frequency  expansion  of 

the  electric  coefficient  b  is  uniform  in  a  and  the  result  for 

n 

perfect  conductivity  can  be  obtained  by  putting  ®  in  the 


expression  for  arbitrary  oenni tti vi ty .  For  the  magnetic  coefficient 
a^ ,  however,  the  expansion  is  not  uniform,  and  for  given  y  ,  ;  f  0 
there  is  a  value  of  ka  Zu^a  below  which  the  solution  is  no  longer 

characteristic  of  a  perfect  conductor.  The  result  for  perfect 
conductivity  cannot  now  be  obtained  from  the  expression  for  a 
general  material,  implying  that,  in  the  present  instance,  there  is 
a  minimum  frequency  at  which  the  facility  can  be  operated.  In 
particular,  a  quasi -magnetostatic  facility  would  be  inappropriate. 

From  (2.53)  the  requirement  is  that 

I kj a|  =  (ka  Zu^o)  >>  1  , 


for  example, 

ka  >  10“  (aZu^a)"^  (2.58) 

If  j  =  10*^  mhos/m  (appropriate  for  silver),  =  1  and  a  =  0.125  m, 
(2.58)  implies. 

ka  >  2.1  X  10"^  , 

corresponding  to  a  minimum  frequency  of  8  kHz.  A  larger  body  would 
permit  a  smaller  minimum  frequency.  If,  in  addition, 

ka  <  0.05 

implying  a  maximum  frequency  of  19  MHz,  the  first  term  in  the  low 
frequency  expansion  is  sufficient,  and  that  term  is  the  one  for  a 
perfect  conductor  (see  2.51).  A  larger  bodv  would  reduce  proportional  1 
the  maximum  allowed  frequency. 
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2.3.2  Effect  of  Penetration  Depth.  It  is  highly  desirable  that  the 
planned  facility  be  able  to  employ  models  which  are  identical  (or  similar) 
to  those  used  in  our  present  facility.  These  models  are,  in  general, 
plastic  with  a  thin  coating  of  silver  paint  to  simulate  a  highly 
conducting  metal.  The  coating  thickness  t  is  typically  0.1  mm  (=  4  mils) 
and  this  should  be  at  least  twice  the  penetration  depth  i  in  the  paint. 

From  (2.53) 


-1/2 


and  the  requirement  therefore  is 

t  >  2  kZu^oj  , 

implying 

k  :  8(Zii^ot2)‘^  . 

If  t  =  0.1  mm,  =  1  and  a  =  10^  mhos/m,  the  minimum  frequency  is 
10  MHz.  This  exceeds  the  minimum  frequency  demanded  in  the  previous 
section,  and  in  conjunction  with  the  maximum  imposed  by  the  require¬ 
ment  that  the  first  term  in  the  low  frequency  expansion  for  a 
0.25  m  model  suffice,  narrows  the  allowed  frequency  range  to 


10  <  f(MHz)  <  19  . 

To  accommodate  a  model  somewhat  larger  than  0.25  m  in  overall 
dimension  (as  well  as  other  considerations),  it  is  suggested  that 
the  facility  be  designed  to  operate  as  close  to  10  MHz  as  possible. 
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2.3.3  Excitation.  If  a  perfectly  conducting  (metallic)  body  is 


exposed  to  an  electromagnetic  field,  the  leading  (magnetostatic)  term 
in  the  low  frequency  expansion  for  the  total  magnetic  field  is  determined 
by  the  leading  term  in  the  low  frequency  expansion  of  the  incident 
magnetic  vector.  In  particular,  neither  the  incident  electric  vector 
nor  the  direction  of  propagation  of  the  incident  field  has  any  effect. 

As  an  example,  in  the  case  of  a  plane  wave  incident  on  a  sphere 
(see  Sec  cion  2.3.1 ) , 
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Of  course,  for  a  more  general  body  will  not  be  simply  proportional 
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The  facility  that  has  been  proposed  consists  of  two  identical 
Helmholtz  coils  of  radius  a  spaced  a  apart.  Between  the  coils  there  is 
a  cylindrical  region  where  the  magnetic  field  is  parallel  to  the  axis 
and  almost  constant  in  amplitude,  with  the  electric  field  circumferential. 
This  is  conceived  to  be  the  test  region  where  the  model  is  placed, 
and  is  in  the  extreme  near  field  of  each  coil.  As  a  result  the  field 
to  which  the  body  is  exposed,  is  quite  different  from  a  plane  wave,  but 
since  the  magnetic  vector  is  constant  in  amplitude  and  direction,  it 
does  simulate  a  plane  wave  as  regards  the  magnetostatic  term.  This  is 
all  that  is  required. 


3.  PHYSICAL  DESIGNS 


As  a  result  of  theoretical  studies  discussed  in  Section  2  it  was 
decided  that  the  facility  would  consist  of  Helmholtz  coils  approximately 
0.5  m  in  radius  spaced  one  radius  apart.  At  zero  frequency  this  would 
provide  a  test  region  of  cylindrical  volume  consisting  of  two  feet  in 
diameter  and  two  feet  in  height  where  the  magnetic  field  is  uniform  to 
within  ±  5  per  cent.  It  has  also  been  shown  that  if  the  current  in  the 
coils  is  constant  in  amplitude  and  phase,  the  frequency  dependent 
variation  of  the  magnetic  field  throughout  the  test  region  will  be  less 
than  one  per  cent  if  f  ^  10  MHz. 

From  the  considerations  of  the  magnetic  field  response  of  typically 
used  models,  the  requirement  on  the  measurement  frequency  is  just  the 
opposite.  A  non-ferrous  metal  does  not  appear  perfectly  conducting  at 
low  frequencies,  and  the  problem  becomes  more  pronounced  for  metallized 
(silver  painted  plastic)  models  used  in  scale  model  measurements.  These 
topics  were  discussed  in  Section  2.3.2  where  it  was  shown  that  for  a 
0.1  mm  (4  mills)  thick  silver  coating  on  a  plastic  model  a  minimum 
frequency  of  aoout  10  MHz  is  required  to  simulate  it  as  a  magnetic 
conductor.  On  the  other  hand,  a  model  made  of  solid  silver  would  allow 
a  minimum  frequency  of  about  8  kHz.  Of  course,  the  use  of  magnetic 
material.;  and  coatings  would  lower  the  minimum  frequency  to  zero  but 
we  have  not  considered  it  at  this  time  due  to  the  difficulties  anticipated 
in  making  or  obtaining  models  of  such  materials.  It  was  therefore  de¬ 
cided  to  construct  a  facility  to  operate  at  10  MHz  but  design  the  feed 


V 


system  for  the  coils  such  that  they  can  be  excited  over  wider  frequency 
range  and,  then,  from  the  results  of  evaluation  measurements  determine 
the  appropriate  operating  range  of  frequencies  for  the  facility.  The 
following  parameters  were  specified  at  the  beginning  of  the  program: 

(i)  geometry  -  coaxial  Helmholtz  coils,  single  turn,  4  feet 
in  diameter  spaced  2  feet  apart 

(ii)  operating  frequency  -10  MHz 

(iii)  excitation  -  at  four  points  on  each  coil 

(iv)  current  in  each  coil  -  as  much  as  practical,  in  the  range 
0.1  to  1 .0  A  (RMS) 

(v)  structure  -  four  4  x  8-foot  plates  butted  together 

With  the  coils  horizontally  placed  as  shown  in  Fig.  3.1  and  using  the  above 
guidelines,  two  facilities  with  different  design  philosophy  were  con¬ 
structed.  One,  using  plywood  support  plates  and  inductive  drivers  for 
the  excitation  of  Helmholtz  coils,  and  the  other,  using  aluminum  support 
plates  and  gap  drivers  in  the  coils.  Henceforth,  the  two  facilities  will 
be  referred  to  as  the  Plywood  and  Aluminum  Facilities,  respectively, 
which  were  constructed  and  tested  also  in  that  order.  The  two  facilities 
are  described  in  Sections  3.1  and  3.2. 

3.1  Plywood  Facility 

3.1.1  The  Structure 

The  facility  was  constructed  by  using  four  4  x  8-foot  3/4-inch 
thick  plywood  sheets  butted  together  as  shown  in  Fig.  3.1  and  held  together 
with  a  frame  made  out  of  2  x  4's.  In  each  sheet  a  semi-circular  hole  of 
radius  \/5~  =  2.24  feet  was  cut  to  form  a  spherical  volume  at  the  center 
that  also  supported  the  coils  and  provided  access  to  the  working  or  test 


Fig.  3.1:  Helmholtz  coil  geometry  with  support 
structure . 
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region.  The  Helmholtz  coils  were  made  from  0.25  inch  (O.D.)  ccpper  tubing 
that  came  in  a  3-foot  diameter  spool  and  hence  was  relatively  easy  to  re¬ 
shape  into  4-foot  diameter  loops  whose  ends  were  then  butted  to  form  a 
short  circuited  loop. 

The  diagram  showing  the  excitation  of  the  facility  is  given  in 
Fig.  3.2.  Originally  the  idea  here  was  to  use  four  equally  spaced 
current  drivers  (such  as,  EG&G  CCDOl  or  CCD02)  on  each  coil,  all  driven 
with  signals  having  equal  amplitude  and  phase.  A  distribution  trans¬ 
former,  was  designed  and  fabricated  to  provide  the  required  eight  uniform 
signals.  The  transformer  and  other  components  associated  with  the  facility 
are  described  in  the  next  section. 

3.1.2  Current  Drivers 

Although  the  CCDOl  or  CCD02  current  drivers  could  have  worked 
satisfactorily,  we  were  not  able  to  obtain  sufficient  numbers  of  them 
and,  therefore,  we  designed  and  built  our  own  drivers. 

An  important  consideration  in  the  design  of  the  drivers  is  the 
load  (Helmholtz  coil)  inductance  which  dictates  the  current  that  can  be 
induced  in  the  coil.  From  computations  and  later  verified  by  measurement 
(Appendix  E)  it  was  found  that  the  inductance  of  a  single  loop  is  about 
5  ;jH,  which  at  10  MHz  has  a  reactance  of  300  ohms  or  75  ohms  per  each 
driver  since  there  are  four  drivers  on  each  coil.  If  the  driver  is  now 
considered  as  a  transformer  with  one  turn  secondary  (a  coil  of  75  ohm 
impedance),  there  will  be  an  impedance  mismatch  at  the  primary  side  when 
referenced  to  50  ohms,  because  more  than  one  turn  is  needed  in  the  primary 
for  the  excitation  of  the  magnetizing  field.  The  torroidal  cores  that 
we  chose  were  T8525  (CMD5005)  from  Ceramic  Magnetics.  Based  on  the 
criterion  that  the  driver  behaves  as  an  ideal  transformer  it  was 


Fig.  3,2:  Excitation  diagram  for  the  Plywood  Facility 


determined  that  six  turns  would  be  reouired  in  the  primary.  Under  this 
condition  and  with  the  secondary  open,  the  measured  impedance  of  the  six- 
turn  primary  was  found  to  be  2000  ohms  at  10  MHz.  (In  retrospect,  it 
now  appears  that  fewer  turns  in  the  primary  may  have  provided  a  better 
overall  performance  of  the  facility,  even  though  the  primary  impedance 
might  have  approached  the  50  ohm  characteristic  impedance  of  the  connec¬ 
ting  coaxial  1 ines . ) 

The  driver  core  with  six-turn  No.  20  copper  wire  winding  was  then 
mounted  in  an  aluminum  housing  shown  in  Fig.  3.3.  Holes  were  drilled  in 
the  housing  to  thread  or  place  the  unit  on  the  coil  which  had  to  be 
broken  (unshorted)  to  do  so.  The  BNC  connector  is  for  connecting  the 
input  line.  Figure  3.4  shows  photographs  of  these  drivers. 

3.1.3  Distribution  Transformer 

The  purpose  of  this  device  was  to  provide  uniform  (equal  amplitude 
and  phase)  signals  to  each  of  the  eight  driver  transformers  mounted  on  the 
two  coils.  Again,  a  torroidal  core  transformer  design  was  used  but  this 
time  with  24-turn  primary  and  eight  3-turn  secondaries.  The  core  selected 
was  Amidon  T-200-2  similar  to  that  used  in  1  KW  amateur  radio  transmitter 
output  circuit.  Even  though  the  choice  of  turns  ratio  was  based  on  con¬ 
siderations  to  provide  uniform  signals  to  each  of  the  eight  outputs,  it 
was  found  that  the  output  voltages  were  not  identical.  The  open  circuit 
voltage  measured  at  7.88  MHz  with  a  10  Meg  probe  on  individual  secondary 
terminals  with  the  other  terminals  terminated  in  100  ohms  varied  from  16.8 
to  1 7  mV  in  amplitude  and  0  to  -2.6  degrees  in  phase. 

Figure  3.5  shows  a  drawing  of  the  aluminum  housing  in  which  the 
transformer  was  placed.  Each  of  the  eight  BNC  outputs  were  connected  with 
matched  lengths  of  RG-58A/U  cable  to  the  respective  current  driver  mounted 


A)  Casing  (  housing  the  toroid  ) 


B)  Toroid  (CMO  3005  T8323  ) 


Fig.  3.3:  Aluminum  housing  for  the  current  drivers. 


(a)  Driver  with  open  cover  to  show  the  core  and  the  winding. 


(b)  Cover  closed.  Note  its  size  relative  to  a  dime. 


Fig.  3.4:  Current  driver  for  Helmholtz  coils.  There  is  a  total 
of  eight  of  these  in  the  facility. 


Fig.  3.5;  Distribution  transformer  housing  for  the 
Plywood  Facility. 
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on  the  coils.  Voltages  were  again  measured,  this  time  with  the  10  Meg  probe 
connected  at  the  output  of  a  given  cable  with  the  other  seven  cables  con¬ 
nected  to  appropriate  current  drivers.  Table  3.1  shows  that  the  amplitude 
of  the  signals  varies  slightly  more  than  the  case  when  measured  at  the 
terminals  of  the  distribution  transformer,  but  the  phase  deviates  by  as 
much  as  9.6  degrees  which  should  be  compared  to  2.6  degrees  obtained  in 
the  previous  case. 

Table  3.1 

Drive  Signal  Voltages 


The  measurements  presented  in  Table  3.1  were  made  using  a  100  mW  signal 
generator  set  at  7.88  MHz.  This  frequency  was  dictated  by  the  optimum 
operation  of  the  faci 1 ity  in  conjunction  with  a  Heathkit  25  watt  amateur 
radio  transmitter  which  we  originally  anticipated  using  as  the  high  power 
source  required  for  making  current  measurements  using  2-3  mm  diameter 
loop  probes.  However,  the  Heathkit  transmitter  was  discarded  after  we 
acquired  a  25  watt  (1  -500  MHz)  power  amplifier. 


3.1.4  Tuning  of  the  Facility 


When  we  started  to  examine  the  current  drivers  for  the  rielmno'tz 
coil  it  became  apparent  that  due  to  impedance  mismatch  in  the  drivers, 
tuning  of  the  facility  would  be  required  to  obtain  efficient  excitation 
of  the  coils.  We  considered  and  tried  tuning  the  source,  the  distribu¬ 
tion  transformer,  the  drivers,  and  even  the  coils  by  introducing  series 
capacitors  in  the  loops.  Most  of  these  approaches  turned  out  to  be  im¬ 
practical  and  we  finally  decided  to  use  sections  of  coaxial  lines  as 
"quarter-wave  matching  devices"  connected  between  the  distribution  trans¬ 
former  and  the  current  drivers.  After  carrying  out  some  experimental  and 
theoretical  (Appendix  D)  studies  the  optimized  length  of  these  cables 
was  found  to  be  18  foot.  When  tuning  the  facility,  either  the  MGL-6A(A)  or 
the  MTL-2A(ER)  sensor  was  used  to  detect  the  mannetic  fields  and  since 

the  two  sensors  provided  sufficient  signal  output  even  with  a  couple 
of  milliwatts  drive,  a  single  generator  (without  an  amplifier)  covering 
1-30  MHz  range  sufficed. 

3 . 2  Aluminum  Facility 

After  it  was  recognized  that  the  Plywood  Facility  with  transformer 
drivers  will  operate  only  at  a  single  frequency  (although  it  could  be 
retuned  by  changing  cable  lengths),  we  also  proceeded  with  the  design 
and  construction  of  the  Aluminum  Facility.  This  facility  uses  gap  exci¬ 
tations  of  coils  instead  of  coupling  transformers  and  metal  plates  instead 
of  plywood  to  ensure  more  uniform  field  distribution  within.  The  metal 
plates  also  allow  routing  of  coaxial  leads  along  the  surface  with  minimum 
interaction  with  the  field. 

The  shape  and  size  of  the  Aluminum  Facility  are  the  same  as  those 
of  the  Plywood  Facility,  and  that  is  where  the  similarities  end.  The 


4  X  8-foot  plates  were  made  of  1/8-inch  thick  aluminum,  the  coils  and 
feed  lines  were  made  of  0.141-inch  semirigid  50-ohm  coax  and  the  connec¬ 
tors  used  were  of  the  SMA-series  type.  The  distribution  transformer 
(North  Hill,  Model  0011)  had  differential  outputs  +\/  and  -V  voltages 
required  to  drive  the  coils. 

Figure  3.6  shows  a  photograph  of  the  facility  with  a  prolate 
spheroid  model  inside.  Dimensions  of  the  facility  are  shown  in  Fig.  3.7. 
The  details  of  implementation  of  the  four  excitation  gaps  in  each  loop 
are  shown  in  Fig.  3.8.  Each  loop  is  fed  by  four  coaxial  lines  originating 
from  the  distribution  transformer.  Four  of  the  lines  (two  from  each  coil) 
are  connected  to  each  side  at  the  output  of  the  transformer  and  to  im¬ 
prove  the  impedance  matching,  39  ohm  (2  watt)  resistors  were  added  in 
series  with  each  line  [8].  The  lines  were  very  carefully  matched  in 
length  to  within  0.125  inches  by  measuring  the  distance  with  a  Time  Domain 
Reflectometer  from  the  input  end  that  connects  at  the  transformer  to  the 
gap  in  the  coil. 
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Fig.  3. 


6:  Aluminum  Facility  shown  with  a  2:1  prolate 
spheroid  (4.5  by  9.0  inches)  mounted  for 
measurement. 
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4.  TESTS  AND  EVALUATION  OF  FACILITIES 


After  construction  of  the  two  facilities,  each  was  evaluated  by 
measuring: 

(a)  the  magnetic  field  to  determine  the  field  intensity  distri¬ 
bution  within  the  coils 

(b)  the  sensitivity,  i.e.,  the  field  intensity  per  given  power 
input 

(c)  frequency  response. 

The  measurements  were  carried  out  by  using  the  experimental  set-up 
shown  in  Fig.  4.1.  An  EG&G  MTL-2A{A)  B-dot  sensor  in  conjunction  with 

DTL-96D  balun  was  used  to  measure  the  magnetic  field.  Since  the  sensor 

2 

has  relatively  large  sensitivity  (A^^  =  0.01  m  ) ,  no  power  amplifier  on 
the  driver  side  nor  pre-amplifier  on  the  sensor  side  was  required.  A 
50-100  mw  CW  signal  generator  was  used  to  drive  the  coils.  Part  of  the 
signal  (from  the  directional  coupler)  was  fed  to  the  vector  voltmeter 
in  Port  A  to  provide  the  reference  signal.  The  test  signal  was  fed 
directly  from  the  sensor  (or  the  probe)  into  Port  B.  Thus,  a  given 
measurement  consisted  of  recording  the  voltage  amplitude  at  Port  B  and 
the  phase  difference  between  the  voltages  at  Ports  B  and  A. 

The  evaluation  measurements  for  both  facilities  were  performed  in 
a  similar  manner.  For  field  probing  the  MTL-2A(A)  sensor  was  scanned 
through  a  cylindrical  volume:  r  =  0,  4,  .  .  .  ,24  inches;  ;  =  0,  45, 

.  .  .  ,  315  degrees;  and  z  =  ±  20  inches  which  resulted  in  measurements 
at  6  X  8  X  11  +11  =  539  points  for  each  test  frequency.  During  these 


measurements  appropriate  care  was  taken  to  position  the  handle  of  the 
sensor  and  the  signal  lead  in  a  radially  outward  direction  so  as  to 
minimize  interaction  with  the  (small)  electric  fields  present  in  the  ; 
directions.  Figure  3.6  shows  the  styrofoam  platform  whose  surface 
contained  the  r  and  ;  coordinates  so  as  to  facilitate  the  placement  of 
the  sensor.  The  entire  platform  was  raised  appropriately  up  and  down  to 
move  the  sensor  in  the  z-di rections . 

4.1  Plywood  Facility  Field  Plots 

Since  the  Plywood  Facility  operates  efficiently  only  at  a  single 
frequency  and  had  to  be  retuned  by  changing  the  power  distribution  cable 
lengths  to  operate  at  another  frequency,  the  field  measurements  were 
carried  out  only  at  7.875  MHz.  As  mentioned  previously,  539  measurements 
were  made  at  a  given  frequency  and  the  data  were  reduced  in  terms  of 
contours  for  r  5S  deviations  of  the  field  from  its  value  at  the  center. 

Figure  4.2  shows  the  results  obtained  in  the  horizontal  planes 
z  =  0  and  z  =  ;  4  inches.  It  can  be  seen  from  Fig.  4.2  that. the  olots 
are  approximately  circular,  and  they  become  larger  as  one  goes  up  or 
down  from  the  center  plane,  since  the  measurement  position  is  then 
approaching  closer  to  the  coils.  Figure  4.3  shows  the  results  obtained 
in  the  horizontal  planes  z  =  t  6  inches,  but  here  the  two  circles  are 
not  of  the  same  size,  and  this  is  due  to  the  fact  that  the  r  5  percent 
boundary  values  vary  rapidly  at  z  =  ±  4  inches  which  can  be  seen  from 
Fig.  4.4  that  gives  the  vertical  plane  results. 

The  vertical  plane  results  in  Fig.  4.4  show  the  uniformity  of 
the  field  intensity  more  explicitly.  The  curves  in  Fig.  4.4  indicate 
that  the  contours  are  not  continuous;  the  +5  percent  contours  are 
near  the  conductors  while  the  -5  percent  contours  are  in  between  the 
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Vertical  plane  results  at  0=0  and  0=90  degrees 
(Plywood  Facility). 


conductors.  As  expected  from  the  circular  plots  obtained  in  the  horizontal 
planes  (Figs.  4.2  and  4.3),  there  is  little  difference  between  the  results 
in  the  0-degree  and  90-degree  planes. 

Overall,  it  was  found  that  the  Plywood  Facility  provides  a  uniform 
field  with  a  maximum  of  ±  5  percent  variations  within  a  spherical  volume 
of  radius  12  inches  and  this  is  consistent  with  the  theoretical  predic¬ 
tions  of  Table  2.1 . 

4.2  Aluminum  Facility  Field  Plots 

Whereas  the  Plywood  Facility  was  probed  at  the  single  frequency  of 
7.875  MHz,  the  field  measurements  with  the  Aluminum  Facility  were  carried 
out  at  7.875,  10.000,  and  15,000  MHz.  As  was  pointed  out  in  the  previous 
sections,  the  horizontal  plane  results  provide  very  little  information 
other  than  the  deviations  from  the  (ideal)  circular  contours;  hence,  we 
show  the  data  for  only  10.000  MHz.  Figures  4.5  and  4.6  show  the  results 
obtained  at  the  horizontal  planes  z  =  0  and  ±  4,  z  *  ±  7  inches  respec¬ 
tively.  The  results  of  Fig.  4.6  show  the  maximum  of  such  deviations, 
and  again,  this  can  be  attributed  to  the  rapid  variations  of  the  fields 
at  this  height. 

Figures  4.7,  4.8,  and  4.9  show  the  vertical  (ip  =  90°)  plane  results 
obtained  at  7.875,  10.000,  and  15.000  MHz,  respectively.  These  results 
indicate  that  the  contour  lines  are  essentially  invariant  with  frequency, 
although  a  closer  inspection  does  show  that  the  :r  5  percent  contour 
tends  to  move  slightly  towards  the  center  with  increasing  frequency. 
Therefore,  we  conclude  that  for  the  frequencies  considered  the  facility 
provides  a  uniform  field  with  a  maximum  of  ±  5  percent  deviation  over  a 
spherical  volume  of  radius  12  inches.  Again,  this  is  consistent  with  the 
design  values  (Table  2.1). 
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On  comparing  the  field  plots  obtained  at  7.875  MHz  with  the  Plywood 
and  Aluminum  Facilities  it  is  found  that  the  field  responses  of  the  two 
are  similar  and,  in  fact,  the  plywood  facility  produces  a  slightly 
larger  working  volume.  However,  this  should  not  be  the  only  criterion 
for  making  a  choice  between  the  two.  A  more  important  criterion  is  the 
frequency  range  over  which  the  facility  can  be  operated  effectively. 

To  determine  the  frequency  range  of  operation,  we  used  the  set-up 
shown  in  Fig.  4.1  with  the  sensor  positioned  at  the  center  of  the  working 
volume,  and  recorded  its  output  voltage  and  drive  voltage  over 
1-30  MHz.  Measurements  were  carried  out  with  the  Plywood  and  Aluminum 
Facilities.  The  frequency  responses  of  the  two  facilities  obtained  from 
the  above  measurements  are  shown  in  Fig.  4.10  which  indicates  that  the 
aluminum  facility  is  superior.  The  field  in  the  Plywood  Facility  is  in 
general  smaller  and  shows  a  resonance  at  12  MHz.  This  frequency  is  a 
function  of  the  coil  and  transformer  inductances  and  feed-line  lengths, 
and  by  varying  the  feed-line  lengths  this  resonant  frequency  can  be 
changed  (see  Appendix  D).  Above  12  MHz  the  Aluminum  Facility  produces 
much  higher  fields  for  the  same  input  drive  and  shows  a  strong  resonance 
at  18  MHz.  Since  at  18  MHz  ka  =  0.23  (k  is  the  wave  number  and  a  is  the 
radius  of  the  coil),  the  peak  response  is  not  due  to  the  loop  resonance 
for  which  ka  =  1  would  be  appropriate.  However,  the  feed  line  lengths 
from  the  transformers  to  the  driver  gaps  in  the  coils  (see  Fig.  3.6) 
are  about  0.25  wavelengths,  and  these  could  be  the  cause  of  the  resonance. 
In  the  measurements  we  have  performed  so  far  the  resonances  have  not  been 
a  problem.  But  if  the  resonance  effects  do  become  critical  in  the  future 
measurements,  these  can  be  reduced  considerably  by  changing  the  matching 


5.  SAMPLE  MEASUREMENTS 


The  evaluation  study  of  the  Plywood  and  Aluminum  Facilities 
indicated  that  the  latter  is  superior,  and  consequently,  the  plywood 
facility  was  dismantled.  The  measurements  presented  in  this  section 
were  all  made  in  the  Aluminum  Facility. 

To  make  the  measurements  of  the  surface  currents  on  a  given  body 
a  sensor  (or  probe)  much  smaller  than  the  MTL-2A(A)  must  be  used.  Designs 
for  miniature  MTL  constructions  were  also  considered  but  not  used.  Then 
we  tried  to  use  our  own  2-3  mm  diameter  loops  that  we  generally  use  in 
the  broadband  (100-4700  MHz)  chamber  measurements  and  found  that  these 
can  indeed  be  used  in  the  present  magnetostatic  measurements.  The  measure 
ments  that  we  made  were  exploratory  in  nature,  and  as  our  measurement 
techniques  improved,  so  did  the  accuracy  of  measurements.  The  models 
used  for  measurement  were  aluminum  spheres,  brass  cylinders,  and  a 
silver  painted  F-106  plastic  model. 

5 . 1  Spheres  and  Cyl i nders 

The  experimental  set-up  used  was  as  shown  in  Fig.  4.1  except  that 
the  MTL-2A(A)  was  replaced  by  a  miniature  (2-3  mm  diameter)  shielded  loop. 
At  a  given  frequency  the  sequence  of  measurements  was  to  record  the  probe 
voltage  with  the  sphere  present  (measuring  total  magnetic  field  H)  and 
without  the  sphere  present  (measuring  incident  magnetic  field  Hq).  The 
ratio  of  the  two  measurements,  i.e.,  H/Hg,  yielded  the  surface  field 
normalized  to  the  incident  field.  Measurements  were  performed  at  8  MHz 


using  loop  probes  in  sizes  from  2.7  to  7.5  mw  diameter  on  3-inch  and 
6-inch  diameter  aluminum  spheres.  The  measured  value  of  H/Hq  ranged 
from  1.35  to  1.88  (from  Mie  theory  the  exact  value  is  1.5).  At  the 
beginning,  the  measurements  were  noisy  and  unrepeatable,  mostly  due  to 
signal  leakage  into  the  signal  cables.  By  carefully  aligning  and  taping 
down  the  leads  in  outward  radial  directions  from  the  magnetostatic 
facility,  we  were  able  to  reduce  most  of  the  coupling  effects;  under 
these  conditions  surface  field  ratios  of  1.45  -  1 .47  were  obtained  with  a 
3  mm  diameter  loop  on  the  6-inch  sphere. 

We  then  proceeded  to  perform  measurements  at  7.1  MHi  using  a  set 
of  brass  cylinders.  A  total  of  11  hollow  cylindrical  shells  (with  open 
ends)  ranging  in  diameter  from  0.125  to  2.125  inches  were  cut  12.0  inches 
in  length  from  thin-walled  brass  tubing.  For  a  measurement  a  given 
cylinder  was  placed  on  its  end  in  the' facility  and  appropriate  readings 
were  taken  with  and  without  the  model.  The  diameter  of  the  loon  probe 
used  was  2.6  mm. 

The  results  of  the  above  measurements  are  shown  in  Fig.  5.1.  It 
should  be  noted  that  in  the  case  of  an  infinitely  long  cylinder  the  sur¬ 
face  field  radio  should  be  2.0.  Our  measurement  results  in  Fig.  5.1  show 
that  the  ratio  is  1.2  for  the  smallest  diameter  and  it  approaches  1.7 
as  the  diameter  increases.  The  1.2  value  obtained  for  small  diameter 
cylinders  is  due  to  the  probe  integration  or  averaging  effects  over  the 
area  of  the  probe.  The  fact  that  the  upper  value  of  1.7  is  different 
from  2.0  may  be  due  to  the  finite  length  of  the  cylinders. 

5.2  Measurements  on  NASA  F-106B  Model 

To  demonstrate  the  usefulness  of  the  facility  for  system  analysis 
studies,  sample  measurements  were  carried  out  with  a  NASA  F-106B  scale 


model  aircraft.  This  was  the  same  1/48  scale  model  used  in  2-100  MHz 
(full  scale)  response  measurements.  With  the  present  facility  we  are 
able  to  provide  data  in  the  20  KHz  to  0.52  MHz  range  using  1  to  25  MHz 
measurement  range. 

These  measurements  were  performed  with  a  25  watt  (1-500  MHz)  power 
amplifier  added  on  the  (coil)  driver  side  and  a  25  dB  low  noise  amplifier 
(1-500  MHz)  introduced  between  the  probe  and  vector  voltmeter.  Also,  the 
signal  generator  was  replaced  by  a  synthesized  source  to  assure  accurate 
repeatability  of  the  drive  signal  frequency  and  the  power  level.  Measure¬ 
ments  were  performed  by  manually  punching  in  the  desired  frequencies  and 
reading  the  probe  voltages  off  the  meter.  This  modified  set-up  did  not 
alleviate  the  problem  of  signal  coupling  through  the  cables  and  connectors 
but,  in  fact,  made  it  worse  due  to  the  added  amplifiers  and  cables.  We 
then  observed  that  the  signal  that  couples  through  the  cables  remains 
essentially  invariant  as  long  as  the  cables  are  undisturbed.  The  polarity 
of  the  probe  signal,  however,  depends  on  the  orientation  of  the  loop. 

Thus,  if  two  measurements  are  made,  one  with  the  loop  oriented  in  a  cer¬ 
tain  way  and  the  other  with  the  loop  rotated  180  degrees,  then  the 
difference  between  the  two  recordings  yields  mainly  the  signal  picked  up 
by  the  loop  and  most  of  the  undesired  signals  are  cancelled  out.  This 
was  then  the  procedure  used  during  the  measurements  with  F-106B  model. 
Figure  5,2  shows  photographs  of  the  model  in  the  facility  set  up  to 
measure  axial  fuselage  current  on  top  of  the  fuselage.  The  top  photograph 
shows  the  model,  the  two  coils  supported  on  plywood  rings,  and  the  al¬ 
uminum  plates.  The  probe  lead  comes  radially  outward,  which  then  is 
taped  to  one  of  the  plates  (not  shown).  The  lower  photograph  in  the 
figure  shows  a  close-up  of  the  model. 


(a)  Facility  and  the  model. 


Figures  5.3  to  5.8  show  the  data  obtained  for  the  following  three 
cases  that  were  also  used  in  previous  anechoic  chamber  measurements  [9]: 

(i )  Top  Inc.,E-parl  fus,  STA:F481R,  J  ;  Figs.  5.3  and  5.4 

d 

(ii)  Top  Inc.,  E-perp  fus,  STA:WL104B,  Ja;  Figs.  5.5  and  5.6 

(iii)  Top  Inc.,  E-perp  fus,  STA:F4816,  Figs.  5.7  and  5.8 

In  each  of  the  above  three  cases,  the  data  are  at  first  presented 
for  the  entire  frequency  range  0-100  MHz  and  there  the  magnetostatic 
result  appears  as  a  dot  at  zero  frequency.  In  the  following  plot  the 
data  are  then  shown  in  the  frequency  range  0  to  3  MHz  where  the  magneto¬ 
static  measurements  are  clearly  seen.  On  looking  at  the  measurements 
presented  it  is  found  that  the  magnetostatic  results  still  have  as  much 
as  1  dB  ripple  in  amplitude  and  15°  ripple  in  phase.  By  developing  better 
measurement  procedures,  for  example,  improving  the  shielding  in  cables 
and  connectors,  and  in  particular,  placing  the  voltmeter  and  signal  ampli¬ 
fier  in  a  shielded  enclosure,  these  ripples  could  be  reduced  to  about 
0.25  dB  and  3  degrees,  respectively. 
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6.  CONCLUSIONS 


The  purpose  of  the  present  program  was  to  develop  an  experimental 
facility  to  measure  the  induced  surface  fields  on  relatively  small 
models  or  objects  at  frequencies  below  the  lowest  frequency  limits 
of  anechoic  chambers.  At  these  low  frequencies  the  magnetic  fields 
(surface  currents)  and  electric  fields  (surface  charges)  are 
independent,  and,  hence,  in  practice  these  can  be  measured  indepen¬ 
dently  in  two  different  types  of  facilities,  e.g.,  magnetostatic  and 
(quasi)  electrostatic. 

Tha  present  program  was  concerned  with  the  development  of  a 
magnetostatic  facility.  To  this  end,  we  conducted  theoretical  studies 
relating  to  the  design  of  the  facility,  and  analyzed  the  responses 
of  perfectly  and  imperfectaly  conducting  canonical  models  place  in 
the  facility. 

The  design,  construction  and  testing  of  two  configurations,  the 
Plywood  Facility  and  the  Aluminum  Facility,  were  carried  out.  The  one 
that  was  selected  and  implemented  was  the  Aluminum  Facility,  which 
met  or  exceded  most  of  the  design  goals.  The  final  dimensions  of  the 
required  Helmholtz  coils  were  based  on  theoretical  computations  to 
obtain  a  cylindrical  working  volume  of  two  feet  in  diameter  and  two 
feet  in  height  having  ±5  percent  maximum  deviation  from  the  field 
at  the  center.  For  the  Aluminum  Facility,  field  mapping  measurments 
were  performed  at  7.875,  10.000  and  15.000  MHz  to  demonstrate  the 
performance  of  the  facility.  Although  working  volume  obtained  from 
these  measurements  met  the  design  goal,  the  data  indicate  that  the 
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volume  is  more  like  spherical  (two  feet  diameter)  rather  than 
cyl indrical . 

No  evaluation  tests  were  performed  above  15  MHz,  but  it  is 
believed  that  the  facility  can  be  used  at  frequencies  up  to  30  MHz, 
perhaps  50  MHz.  In  the  original  guidelines  of  the  program  the  intent 
was  to  operate  the  facility  at  10  MHz,  but  by  using  broadband  components 
we  have  demonstrated  its  capability  through  the  range  of  frequencies 
of  at  least  1-25  MHz  by  making  sample  measurements  on  F-106B  model 
aircraft. 

To  provide  a  complete  set  of  low  frequency  electromagnetic 
response  data  for  a  given  model,  it  is  suggested  that  the  magnetostatic 
facility  be'  complemented  by  a  (quasi)  electrostatic  facility  capable 
of  measuring  induced  surface  charges  on  the  model.  Theoretical,  yet 
practical,  designs  have  been  considered  for  such  a  facility  which  would 
be  similar  in  size  to  the  magnetostatic  one,  but  with  the  Helmholtz 
coils  replaced  by  conductive  spherical  shells  of  comparable  dimension. 
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APPENDIX  A 


ADDITIONAL  DETAILS  FOR  SECTION  2.1 

A.l  Derivation  of  Eqs.  (2.9)  and  (2.10) 

In  this  Appendix  we  show  the  derivation  of  the  approximate  expressions 
(2.9)  and  (2.10)  for  and  from  the  corresponding  exact  expressions 

(2.7)  and  (2.8)  which  are  rewritten  as  follows; 
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(A. 3) 


and  K  and  E  are  defined  in  (2.4)  and  (2.5).  Note  that  (A. 3)  indicates 
k  *  0  as  r  -  0.  for  all  z  since  a  0,  and  k"  <  1  for  r  ^  0  and  for  all  z. 
We  shall  make  use  of  the  following  known  expansions  of  K  and  E 

(see  [6] ' : 
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A. 2  Derivation  of  Eg.  (2.9) 

To  obtain  an  expression  for  near  r  0  we  rewrite  (A.l) 
as : 
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Expanding  (1  -  k2/2)/(l  -  k2)  as  a  power  series  in  k2,  and  using 
(A. 4)  and  (A. 5)  it  can  be  shown  that 
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Now  substituting  (A. 7)  in  (A. 6)  we  obtain 
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A, 3  Derivation  of  Eg.  (2.10) 

The  expression  for  near  r  -►  0  is  obtained  from  the  Taylor's 
series  expansion  (A. 2)  expressed  as: 
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where  and  its  derivatives  on  the  right-hand  side  are  obtained 
from  (A. 2)  as  discussed  below. 
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Using  (A. 4)  and  (A. 5)  it  is  simple  to  show  from  (A. 2)  that 
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To  obtain  the  remaining  terms  on  the  right-hand  side  of  (A. 9)  we 
need  the  following  relationship  obtained  from  (A. 3): 
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and  the  following  derivatives  of  K  and  E  ^unctions  at  r  *  0  obtained 
from  (A. 4)  and  (A. 5)  in  conjunction  with  (A.ll); 
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After  differentiation  of  (A. 3)  and  using  the  appropriate  relations 
in  (A. 12)  it  can  be  shown  after  considerable  algebraic  manipulation 


that: 
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How,  introducing  (A. 13)  into  (A. 9),  we  obtain 
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wHicn  is  the  required  equation  (2.10). 


The  following  is  a  listing  of  a  program  (entitled  SQ81:  MELMITER) 
for  computing  the  fields  within  the  Helmholtz  coil  geometry  and  based  on 
Eqs.  {Z.lt)  ana  (2./d7). 


1 

IMPLICIT  REAL48  <A-H.J-Z) 

-) 

INTEGER44  I,MtN 

3 

FZ< I )=Z+I-1 .5 

4 

K(  I)=DSaRT(4#R/<  <RH  )**2+FZ<  I  )**2>  > 

5 

FREF<  I  )=5.#IiSQRr<5.D0)/16./PI/DSQRT(  ( R+ 1  )  1i*2+FZ (  I  )  **2) 

6 

ELINTE<  I  )*<  I-l)*ELIF-E2+<2-I  )*ELIPE1 

ELINTKl I )=( [-1 )*ELIPK2+(2-I XELIPKl 

9 

HR< I )=PREF( I )*F2( C )/R#<-EL INIM I >+( l+RUR+FZC I ) *42 ) / 

9 

2  ( (R-1 >**2+FZ< I )4*2)#ELINTE< I >  ) 

10 

HZ( I)=PREF( I )*<ELINTK( I )  +  ( 1-R*R-FZ( I  )**2)/ 

11 

2  ( <R-1 )»*2+FZ( I )4*2)*ELINTE(  I  )  ) 

12 

PI=3. 141592653589793238D0 

13 

1 

URrrE<6. 11 ) 

14 

11 

FORMAT('  enter:  R1 . Z1 f R2 » Z2 »N . RI 1 » Z I 1 , RI 2 » ZI 2 f M . ' ) 

13 

REA0<5. 12)R1 .Z1 .R2.Z2fN.RIl.ZIl »RI2»ZI2fM 

16 

12 

F0RMAT<2<4F10.5. 15) ) 

17 

klRITE<6>13) 

18 

13 

FORMAT! .9Xt 'R/A' » 12X» 'Z/A' » 13Xf 'HR/HO' » 16X. 'HZ/HO' / ) 

19 

IF(N.LE. 1 )N*1 

20 

IF(M.LE.1 )M»1 

21 

ZINC»0 

22 

ZlINC-0 

23 

RINC*0 

24 

RI1NC=0 

25 

IFIN.EO. 1 'G0T021 

26 

ZINC  =  <Z2-Z1 )/(N-l  ) 

27 

RINC=<R2-R1 )/<N-  1  ) 

29 

21 

IFIM.EO. 1 >001022 

29 

ZIINC=<ZI2-ZIl>/<M-l) 

30 

RIINC  =  (RI2-RI1 )/<M-T  ) 

31 

-n 

Z1=Z1+ZI1 

32 

R1=R1+RI1 

33 

DO  102  12=1. M 

34 

Z  =  Z1 

35 

R=R1 

36 

DO  101  11  =  1. N 

37 

IF  (R.LE. 1 .D-20)R=1 .D-20 

38 

ELIPEl=DeLIEl<R<l)) 

39 

ELIPE2=DELIE1<K(2)  ) 

40 

ELIPK1  =  DELIK1<K<1)  ) 

41 

ELIPK2=0ELIM (K(2) ) 

42 

HRTOT=HR( 1 )+HR<2> 

43 

HZTOT-HZ< 1 )+HZ(2) 

44 

WRITE<6. 14)R.Z.HRTnT.HZT0T 

45 

14 

F0RMAT(2F15.6.2F21 .10) 

46 

Z=Z+ZINC 

47 

101 

R=R+RINC 

48 

Z1=Z1+ZIINC 

49 

102 

R1=R1+RI INC 

50 

URITE<6.15) 

51 

15 

FORMAT!// ^//> 

32 

STOP 

33 

END 
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APPEXDIX  B 


MAGNETOSTATIC  FIELDS  NEAR  A  CONDUCTING  BODY 


The  formulation  and  sample  computations  of  the  magnetostatic 
fields  on  and  near  a  perfectly  conducting  sphere  and  a  prolate  spheroid 
are  presented.  In  the  static  and  quasi-static  regimes  a  prolate  spheroid 
could  be  used  to  approximate,  for  example,  an  aircraft  fuselage.  Giri 
and  Sands  [B.l]  have  used  somewhat  simplified  expressions  to  predict 
the  interaction  errors  for  the  frontal  region  of  the  NASA  F-106.  Because 
the  errors  they  predict  depend  on  the  particular  model  chosen  to  represent 
the  aircraft,  plus  the  fact  that  the  scale  model  measurements  show  that 
most  of  the  interaction  errors  occur  in  aircraft  resonance  regions  where 
quasi-static  representations  are  invalid,  we  have  chosen  not  to  pursue 
the  generation  of  data  that  would  have  questionable  relevance  to  the 
problem.  Instead,  we  present  here  rather  general  magnetostatic  analyses 
for  the  sphere  and  the  prolate  spheroid. 


B. 1  Sphere 

For  a  sphere  of  radius  a  illuminated  by  the  plane  wave 


X  e 


-ikz 


the  exact  (dynamic)  solution  is  available  in  the  form  of  a  Mie  series. 
From  Bowman  et  al  [B.2],  the  total  (incident  plus  scattered)  magnetic 
field  components  are 
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where  a  time  factor  6’’“^  has  been  assumed  and  suppressed.  The  notation 


is  defined  in  the  above  reference. 

To  the  leading  order  in  k  only  the  terms  corresponding  to  n  »  1 
contribute.  Since 
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the  magnetostatic  field  components  are 


and  thus 
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-Y  cos 
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^B.l 


These  are,  of  course,  independent  of  the  direction  of  incidence  of  the 
original  plane  wave,  and  the  normalized  tangential  components  are  equal 
to  3/2  on  t  e  surface  whereas  the  radial  component  is  zero  there. 


B.  2  Prolate  Spheroid 

It  is  now  simpler  to  obtain  the  magnetostatic  fields  directly 
rather  than  as  the  limit  of  a  dynamic  solution. 

In  terms  of  the  prolate  spheroidal  coordinates  t,",;  where 
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with  1  ^  •  'o,  -1  1  1  and  0  ^  ;  <  2-,  the  spheroid  is  defined  as 

the  surface  ;  =  As  evident  from  these  relations,  the  z  axis  is  the 


axis  of  rotation  and  d  is  the  semi-interfocal  distance.  The  length  to 
width  ratio  of  the  spheroid  is 
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J-1/2 


and  ranges  from  «  for  a  long  thin  spheroid  or  needle  to  1  for  a  sphere. 
The  corresponding  values  of  ;;  are  1  and  «  respectively. 


(i)  If 


the  incident  magnetostatic  potential  is 
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We  seek  a  scattered  magnetostatic  field  =  -74/^  where  is  an 
exterior  potential  satisfying  the  boundary  condition 


on  the  surface.  Since 
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we  are  led  to  choose 
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and  application  of  the  boundary  condition  then  gives 
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In  terms  of  the  shape  factor  q  [B.3]  where  q  varies  monotonically 
from  2/3  for  a  sphere  to  1  for  a  needle. 
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The  same  results  are  obtained  if  h’  =  y, 
(ii)  If 
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with 
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the  form  assumed  for  the  scattered  magnetostatic  potential  is 
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V, 

V, 
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and  application  of  the  boundary  condition  then  gives 
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In  terms  of  the  same  shape  factor  q 
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B. 3  Numerical  Results 

For  a  sphere  any  tangential  component  of  the  total  magnetostatic 
field  when  normalized  to  the  corresponding  component  of  the  incident 
field  decreases  from  a  value  1.5  at  the  surface  to  1  at  large  distances, 
whereas  the  radial  component  when  similarly  normalized  increases  from  0 
at  the  surface  to  1  at  large  distances.  If 


r  =  a(l  +  y) 

% 

so  that  Y  is  the  distance  from  the  surface  expressed  as  a  fraction  of 
the  radius,  (B.l)  and  (B.k!)  give 


nor 


=  ^  =  1  -  (1  +  y) 


-3 


‘  tan  = 
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and  these  are  illustrated  in  Fig.  B.l.  We  observed  that  the  percentage 
reduction  in  a  tangential  component  below  its  surface  value  of  1.5  is 
approximately  the  same  as  y,  i.e.  the  distance  r-  a  expressed  as  a 
percentage  of  the  radius. 

The  results  for  a  spheroid  are  more  difficult  to  express  comoactly, 
and  since  the  normalized  field  components  depend  on  the  orientation  of 
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the  spheroid  relative  to  the  incident  magnetic  field,  it  is  convenient 
to  consider  separately  the  cases  when  H  *2=0  and  H  •  z  =  1.  These 

will  oe  referred  to  as  the  perpendicular  (or  transverse)  and  parallel 
or'entations  respectively.  ,  rom  (B.3)  -  (B.6)  we  then  have 
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and  these  are  plotted  as  functions  of  the  length  to  width  na' 
spheroid,  1  ^  i/vt  ^  10,  in  Fig.  B.2. 

To  show  the  variation  of  the  normalized  ^ielos  as  a  ■ 
the  distance  away  from  the  surface,  we  nave  p.a-  -f--  '  • 
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spheroids  having  z/w  =  2  and  5,  for  which  =  1.15470  and  1.02062 
respectively.  In  the  plane  2=0  the  distance  from  the  surface  expressed 
as  a  fraction  of  the  semi  minor  axis  is 


.^nd  for  both  spheroids  we  have  computed  r  and  r.  as  functions  of 

nor  tan 

Y,  0  <^Y  1,0.1.  The  results  are  shown  in  Figs.  B.3  and  B.4,  and  as  in 
the  case  of  the  sphere,  the  results  for  the  tangential  components  are 
given  as  percentage  reductions  below  the  surface  values.  We  observe 
that  the  increase  in  away  from  the  surface  is  almost  the  same  for 
both  orientations  of  the  spheroid,  and  the  rate  decreases  slightly  with 
increasing  a/w.  For  the  tangential  components,  the  percent  reduction 
is  almost  independent  of  «-/w  for  the  perpendicular  orientation,  whereas 
for  the  parallel  orientation  the  reduction  decreases  rapidly  as  £/w 
increases.  Indeed,  the  decrease  is  roughly  proportional  to  (£/w)2. 
Finally  we  remark  that  y  is  only  a  measure  of  the  fractional  distance 
from  the  surface  in  the  meridional  plane  z  =  0.  As  we  move  round  the 
spheroid  towards  the  tips,  y  translates  into  a  smaller  physical  distance 
from  the  surface,  and  at  a  tip,  the  distancr  is  smaller  by  a  factor  z/w 
(approx.)  and  by  a  factor  (ii/w)^  v;hen  expressed  as  a  fraction  of  the 
same  major  axis. 
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Fig.  B.3:  ( — )  and  percent  reduction  in  ( - )  for  a  2:1 

prolate  spheroid  as  functions  of  distance  away  from  the 
surface. 
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APPEND  IK  C 

LOU  FREQUENCY  SURFACE  FIELD  BEHAVIOR 

The  main  reason  for  developing  a  quasi-static  measurement 
facility  is  to  better  define  the  transfer  functions  at  frequencies 
below  those  for  which  the  present  surface  field  measurement  facility 
is  effective.  As  the  model  frequency  is  decreased  below  (about) 

150  MHz  the  performance  of  the  anechoic  chamber  deteriorates  to  such 
an  extent  that,  in  combination  with  the  approaching  cut-off  of  the 
illuminating  antenna,  the  measured  data  for  the  surface  fields 
(currents  and  charges)  can  be  significantly  in  error;  and  to  achieve 
the  same  full  scale  frequencies  by  using  a  smaller  model  would 
require  working  with  a  model  that  was  too  small  to  be  conveniently 
handled,  and  with  a  probe  that  was  no  longer  small  compared  to  the 
dimensions  of  the  structure  on  which  the  fields  were  measured. 

At  frequencies  for  which  the  model  is  small  compared  to  the 
incident  wavelength,  any  surface  field  quantity  V  can  be  expanded  in 
the  form 

V  *  A  +  iu,B  +  uj2c  +  0(oj3)  (C.l ) 

where  w  is  the  circular  frequency  corresponding  to  a  time  convention 
and  A,B,C  are  real  functions  of  position  or.  the  model.  If 
A  0,  then 
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for  sufficiently  small,  and  in  a  neighborhood  of  tu  =  0  the  resulting 
curve  is  parabolic:  concave  up  if  A  and  C  +  0.5B^/A  have  the  same 
sign,  and  concave  down  if  the  signs  are  opposite.  If  A  =  0, 


IVI  -  1B|  (C 

and  the  curve  is  linear,  starting  from  zero.  These  conclusions  are 
independent  of  the  choice  of  phase  origin. 

There  are  a  few  simple  bodies  for  which  A,  B,  and  C  can  be 
obtained  from  a  knowledge  of  the  exact  solution,  and  others  for  which 
they  can  be  found  by  application  of  low  frequency  techniques.  It  is 
of  interest  to  examine  the  results  for  ;ome  examples  of  the  former 
class. 

C.l  Acoustically  Hard  Sphere 
For  the  plane  wave 


.,inc  ikz 
V  =  e 


(C 


incident  on  an  acoustically  hard  sphere  of  radius  a,  center  at  the 
origin  of  coordinates,  the  total  field  at  the  surface  is  [C.l,  p.  374] 
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where  p  =  ka.  Since  [C.3] 
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we  have 
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where  n  =  cos  9,  varying  from  1  at  the  front  of  the  sphere  to  -1  at 
the  back.  Thus 
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and  the  amplitude  curve  is  concave  down  for  29®  <  9  <  151®,  but 
concave  up  otherwise.  For  a  vanishingly  small  sphere  the  incident 
field  alone  satisfies  the  boundary  condition,  and  it  is  therefore 
logical  that  V(0)  =  1. 


C.2  Acoustically  Hard  Prolate 

To  obtain  some  feeling  for  the  effect  of  the  body's  shape, 
consider  now  a  prolate  spheroid  illuminated  by  the  plane  wave  (C.4)  at 
nose-on  incidence.  From  [C.4]  the  total  field  at  the  surface 
is 
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and  since  [C.2,  p.  32] 
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where  c  =  ka,  a  being  half  the  interfocal  distance. 

Using  the  expansions  in  the  Appendix  of  [C.4]  it  can  be  shown 

that 
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By  inserting  the  expressions  for  the  angular  Legendre  polynomials 
in  terms  of  n  =  cos  0  (see,  for  example,  [C.5,  p.  608])  and  using 
the  facts  that 
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we  finally  obtain 
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As  expected,  V  =  1  when  c  =  0,  and  the  expression  reduces  to  (C.6)  when 
C  ->■  «,  c  -►  0  in  such  a  way  that  c^  tends  to  the  finite  value  ka. 

Near  to  the  ends  of  the  spheroid  when  n  =  ±1 ,  the  amplitude 
curve  is  always  concave  up,  but  as  in  the  case  of  the  sphere  there  is 
a  range  of  n  (and  hence  e)  where  the  curve  changes  to  concave  down, 
and  such  a  range  exists  for  any  spheroid.  Thus,  for  S  =  2  corresponding 
to  a  1.15:1  spheroid,  the  curve  is  concave  down  for  26.5°  <  s  <  153.5°. 
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C.3  Perfectly  Conducting  Sphere 

With  a  perfectly  conducting  body  the  surface  field  quantities 
that  are  measured  are  the  normal  component  of  the  electric  field  (or 
charge)  and  the  tangential  components  of  the  magnetic  field  (or  current). 
It  would  be  desirable  to  generate  the  low  frequency  expansions  of 
these  quantities  for  a  perfectly  conducting  spheroid  and,  possibly, 
for  other  simple  shapes  as  well,  but  since  this  is  a  non-trivial  task, 
we  shall,  for  the  moment  at  least,  confine  attention  to  a  sphere. 

For  a  sphere  of  radius  a  illuminated  by  the  incident  plane 
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tne  non-zero  components  of  the  (total)  surface  field  are: 
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The  amplitude  curve  for  is  always  concave  up.  That  for 
Hg  is  concave  up  throughout  the  illuminated  region,  i .e. ,  0  <_  0  <  90' 
and  the  curve  for  H  is  concave  up  except  for  51®  <  0  <  129®. 


C  . 4  Some  Implications 

It  is  anticipated  that  the  quasi-static  facility  will  provide 
data  over  a  frequency  range  comparable  to  or  exceeding  that  for  which 
the  above  low  frequency  expansions  are  valid,  and  will,  in  effect, 
extend  the  frequency  coverage  of  our  present  data  down  to  zero.  When 
that  facility  is  in  operation,  it  will  be  of  interest  to  see  if  the 
low  frequency  data  can  be  adequately  approximated  using  the  theoretical 
values  for  a  simplified  structure- model ing  the  surface  where  the 
measurement  is  made. 

The  main  purpose  for  the  facility  is  to  expand  the  frequency 
coverage  of  the  measured  transfer  function  data  used  for  EMP 
extrapolation,  but  it  is  possible  that  low  frequency  information  could 
also  help  to  determine  the  lowest  (zeroth  order)  SEM  pole.  To 
explore  this,  consider  some  transfer  .function  F(s)  where  s  =  iu.  At 
low  frequencies,  F(s)  can  be  expanded  as 


F(s)  =  F(0)  +  sF'(O)  +  F"(0)  +  .... 


where  the  coefficients  in  the  series  are  real  functions  of  position, 
orientation,  etc.,  as  well  as  of  the  body  itself.  The  SEM  represen¬ 
tation  of  F(s)  is 


F(s)  .  2 


where  the  asterisk  denotes  the  complex  conjugate  and  if  the  poles  s^ 
are  ordered  in  increasing  magnitude,  we  have,  for  |s|  <<  |Sq|, 
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APPENDIX  D 


ANALYSES  OF  EIGHT-PORT  NETWORK 

These  analyses  were  performed  to  predict  and  optimize  the  performance 
of  the  Plywood  Facility  (Section  3.1). 

D. 1  Theoretical  Analysis 

The  analysis  of  the  network  is  based  upon  the  following  assumptions 
about  the  driver  and  coil  circuit  (see  Fig.  D.l): 

A)  Currents  I,  I'  are  the  same  in  all  eight  secondary  distribution 
loops. 

B)  The  mutual  inductance  on  the  primary  distribution  transformer 

depends  only  on  the  relative  positioning  of  each  winding  on  the  toroid. 
Thus  M  =  M-j  +  M2  +  .  .  .  +  My  =  M2  +  •  •  •  +  Mg,  etc. 

C)  All  real  resistances  are  assumed  very  small  and  neglected. 
Definition  of  Symbols; 

M^  =  mutual  inductance  between  primary  and  secondary  windings  on 
primary  distribution  transformer. 

Mg  =  mutual  inductance  between  and  Lg 
M^  =  mutual  inductance  between  helmholtz  coils 
M  =  mutual  inductance  between  secondary  windings  on  primary 
distribution  transformer 
Ly  =  self  inductance  of  primary 
*  inductance  of  secondary 
Lp  =  inductance  of  primary  feed  transformer 


L,  =  i nductance  of  He! r.ihol  tz  coil 


M„  =  mutual  Inductance  of  Helmholtz  coil 


Vj  =  input  voltage 


I,  =  current  in  one  Helmholtz  coil 


We  wish  to  derive  a  transfer  function  which  relates  the  input 


voltage  Vy  to  the  output  current  The  dependent  variables  are  the 


coaxial  cable  length  and  the  frequency  of  operation.  The  analysis  is  done 


in  the  frequency  domain  with  assumed  exp(jcot)  time  convention. 


We  begin  by  analyzing  one  of  our  eight  secondary  distribution  sub¬ 


systems,  as  shown  in  Fig.  D.2. 


(x) 


1  -  r  r.e 

g  I 


e-Ji3x  ^  g-j25^gjex 
I 


(D.l  ) 


where 


\ 


*  ^0 


■  'g  -  ^0  •  ^0 


The  other  case  we  have  is  that  of  Fi gure  D.  3,  essential  1  y  the  mirror  image. 


Replacing  x  with  {i  -  x)  we  then  have 


V2(x)  = 


l  7  +  Z 

'■  ^0  L 


1  -  r  r.e 
g  i 


-j22Z 


e-j3i  gjBx  ^  .  ^-j3x 

g 


(D.2) 


The  resultant  superposition  of  these  two  solutions  (Fig.  D4)  is  our  desired 


resul t : 


\i  ^r\® 


-jBi  JBx  -jSx  , 


--jBx  ,  „-j23i_jBx 


The  system  of  equations  describing  our  whole  circuit  is  then; 

1)  Vg  -  V(^)  +  sL^I  =  0 

2)  V(0)  +  V|_  +  sLgl'  =  0 

3)  s4Mor  -  sM,I  ,  +  sL,  I  ,  =  0 

B  2,  X,  I  i. 

4)  sL^I^  +  sBM^I  -  =  0 


(D.4) 


The  voltages,  which  are  a  function  of  the  coaxial  cable  length,  can  be 
conveniently  written  as: 


where 


V(0)  =  AI^  +  BI  +  Cl 


(D.5) 


^"a^o  ^  -q) 

”o  *  -  Vo'"'''*' 


B  =  -  (1  *  r„) 

<^0  <'  - 


^"b^o  ^ 

<^0  *  -  ^g‘o=''^®''> 


=  DI-p  +  El  +  FIj^ 


(0.6) 


where 


sM^z„  (1  t 

1^0  "  Zg)  (1  -  Tg  r„e-J28‘) 
-  sMZg  (1  t  rje'J^®'') 

(Zq  "  ^g>  <’  - 


F  ,  sM„2g  e-Jg"  (1  .  r„) 

<^0  *  \>  <’  -  ■g''o^'^“''> 


The  simultaneous  solution  of  these  equations  yields  our  desired  transfer 


function; 


^  {sBM«  -  BD  -  sAL,  +  sAM  +  AE} 


+  sMq  +  C  +  F(B  -  fl 

-  «  I  T7 


D.2  Results 

This  transfer  function  was  programmed  and  plots  were  made  for  given 
frequencies  and  coaxial  cable  lengths.  Values  for  the  assorted  self  and 
mutual  inductances  were  either  measured  or  calculated  (where  measurements 
were  deemed  impractical).  They  are: 

=  4.5  X  10'^  H 

=  7  X  10"^  H 

M  =  4.2  X  10'^  H 

Lg  =  43  X  10'®  H 


=  2.8  X  10"  H  (calculated) 

L,^  =  4.6  X  10"^  H  (calculated) 

L.^  =  6.6  X  10'^  H 
Mg  =  5.8  X  10"^  H 

The  computer  program  listing  follows,  along  with  a  computed  response 
(Fig.  0.5)  for  a  14-foot  long  feed  cable.  The  measured  data  are  given  in 
Figs.  0.6  and  0.7  for  18  ft.  and  14.3  ft.  feed  cables,  respectively.  For 
unexplained  reasons  the  18-foot  cable  length  data  matches  the  theory  better 
than  the  14-foot  cable  length  data,  and  at  this  time  we  do  not  have  an 
explanation  for  the  discrepancy. 


Listing  of  the  Program 


1 

real  MIL»LEN.LHfLLrLA»LTfL 

n 

REAL  MArM.HL»M£i 

3 

i; OMRLEX  G  AMO  ,  GAMGO » X  ,  PEN  f  OMEG .  A  r  P  .  C  -  P  >  E  r  F 

A 

COMPLEX  T1 ,T2.T3.T4tT5f IL 

tr 

COMPLEX  Z.ZZ.CMPLX 

6 

COMPLEX  NUM.PENOM 

7 

LA-4.5E-7 

8 

MA=7.0E-7 

9 

M=4.2E-7 

10 

LB=43.E-6 

11 

ML=2.8E-7 

12 

LL=4.6E-6 

13 

LT=6.6E-6 

14 

M£l»5.8E-6 

15 

PI=3. 14159 

16 

WRITE(6.50) 

17 

READ  (5.20)  FREQ.LEN.VT 

18 

Y=2«PI»FREQ 

19 

0MEG=CMPLX(0.0» Y) 

20 

BETA=2*PI*FREQ/(0.666*3.0E8) 

21 

L=LEN*0.304 

X=QMEG*LA 

23 

GAMG=(X-50. )/(XP50. ) 

24 

X=OMEG*LB 

25 

GAM0=(X-50. )/<X+50. ) 

26 

U=C0S<BETA*L) 

27 

0=SIN(-BETAfL) 

28 

UU=C0S(2*BETA«L) 

29 

00=SIN<-2»PETA*L) 

30 

Z=CMPLX<U.y) 

31 

zz=CMPLX(uu»yv) 

32 

PEN»1 .0-GAMG*GAM0*ZZ 

33 

X=50.*c 1 .0+GAMO)«Z/(5O.+0MEG*LA) 

34 

A=0MEG*MA*X/DEN 

35 

B=OMEG»M*X/PEN 

36 

X=50 . *  < 1 , 0  ^GAMG*ZZ ) / ( 50 . +OMEG*LB  > 

37 

C=OMEG»MB*X/DEN 

33 

X=50.*( 1 .0+ZZ*GAM0)/(50.+0MEG*LA) 

39 

D=0MEG»MA*X/DEN 

40 

E=OMEO*M*X/PEN 

41 

X  =  50  .*Z»( 1 .0+GAMG)/(50. tOMEG*LD) 

42 

F=OMEO*MB*X/DEN 

43 

T1=0MEG»B*MA 

44 

T2=B«P 

45 

T3=0MEG*A«LA 

46 

•^4  =  0MEG*A*M 

47 

T5=A«E 

48 

NUM=VT»<T1-T2-T3+T4+T5)/(0MEG*LT ) 

49 

T 1 =OMEG*LB* ( ML-LL ) / ( 4 . 0*MP ) +0MEG4MP+C 

50 

T2=< -0MEG*8.*MA*MA/LT+8.*MA*P/LT-0MEG*M+0MEG*LA-E ) 

51 

T3=(B-8.0»A*MA/LT)*F 

52 

DEN0M=T1*T2+T3 

53 

54 

IL=NUM/PENOM 

55 

MIL=CABS< IL)*1000 

56 

ANGIL=ATAN2< AIMAG< IL) rREAL( IL) ) 4180/PI 

57 

PRINT  100.MIL 

58 

PRINT  200.ANGIL 

59 

PRINT  300. LEN 

60 

PRINT  400. FREQ 

61 

PRINT  500. VT 

62 

20 

FORMAT (3G9. 1 ) 

63 

50 

FORMAT('  ENTER  FREQUENCY  IN  HZ. LENGTH  IN  FT.L'OLTAGE 

64 

100 

FOPMA'i  <  •  COIL  CURRENT=  '  .F7. 2.  '  MILLIAMPS') 

65 

200 

FORMAT('  PHAGE=' .F6. 1 . '  DEGREES') 

66 

300 

FORMAT ('  COAXIAL  CABLE  LENGTH= ' . F7 . 2 . '  FT') 

67 

400 

FORMAT ('  FREQUENCY^' ,F11 . 1 , '  Z') 

68 

500 

FORMAT ('  INPUT  OOLTAGE ^ ' .F6 . 1 . '  OOLTS') 

69 

STOP 

70 

END 
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APPENDIX  E 


INDUCTANCE  OF  A  LOOP 

A  single  turn  loop  was  constructed  to  confirm  the  theoretical 
value  of  inductance  calculated  from  the  formula  [E.l] 

L  =  [7.353  log^Q  {^)  -  6.386]  (E.l) 

where  a  is  the  loop  radius  and  d  is  the  wire  diameter.  With  the  values 
a  =  24  in.  and  d  =  0.125  in.  we  have  L  =  4.62  uH. 

To  test  this  result,  a  four-foot  diameter  loop  was  constructed 
using  RG-58/U  coax  mounted  on  a  styrofoam  sheet  as  shown  below. 


Measurements  were  made  using  a  Boonton  Radio  Corp.  Q-meter,  Type  260-A, 
and  from  the  readings  we  determined  the  Q  and  the  resonant  capacitance  C 
in  the  frequency  range  3.5  to  10.0  MHz.  Inductance  was  found  from  the 


equation  =  1/vAX,  and  the  series  and  parallel  resistances  were  calcu¬ 
lated  for  the  equivalent  parallel  and  series  circuits  shown  in  Fig.  E2. 


Rp  =  (ijLQ  =  cjL/Q 


Fig.  E.2:  Equivalent  circuits. 
Table  El  below  lists  the  computed  equivalent  components. 


Table  El.  Measured  Loop  Parameters. 


f(MHz) 

Q 

C  (pF) 

L  (yH) 

Rg-  (rO 

Rp(kQ) 

3.5 

140 

424 

4.88 

0.767 

15.0 

5.0 

108 

198 

5.12 

1 .49 

17.4 

6.5 

98 

no 

5.45 

2.27 

21  .8 

7.9 

60 

68 

5.97 

4.94 

17.8 

10.0 

18 

33 

7.68 

26.8 

8.69 

Note  that  as  the  frequency  decreases,  L  approaches  the  4.62  uH  value 
computed  from  (E.l). 
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